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Abstract 

Based on a detailed modeling of the Fermi surface topology of MgB2 we calculated the anisotropy of the upper 
critical field B C 2 within the two gap model. The cr-band is modeled as a distorted cylinder and the 7r-band as a 
half-torus, with parameters determined from bandstructure calculations. Our results show that the unusual strong 
temperature dependence of the B C 2 anisotropy, that has been observed recently, can be understood due to the small 
c-axis dispersion of the cylindrical Fermi surface sheets and the small interband pairing interaction as obtained from 
bandstructure calculations. We calculate the magnetic field dependence of the density of states within the vortex 
state for field in c-axis direction and compare with recent measurements of the specific heat on MgB2 single crystals. 
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It now appears widely accepted that the supercon- 
ducting state of MgB2 has to be described by two differ- 
ent s-wave gaps [1] . This also implies unusual behavior 
in its vortex state, which has been observed recently. 
The anisotropy of the upper critical field B C 2 displays 
an anomalous strong temperature dependence. An up- 
ward curvature of B C 2 close to T c is observed, when the 
field is directed in ab-plane direction [2,3]. The specific 
heat and thermal conductivity show a rapid increase as 
a function of magnetic field [4,5] and the size of the vor- 
tex core and local density of states show unusual behav- 
ior as well [6] . While it seems reasonable that these ob- 
servations are related to the two-gap nature of MgB2 , 
a detailed understanding of all of these phenomena 
within a single model is still lacking. We have shown 
recently that the unusual behavior of the upper critical 
field can be understood, if the topology of the Fermi 
surface and the weak interband pairing interaction is 
taken into account [7]. Within our model the cr-band is 
described as a distorted cylinder and the 7r-band as a 
half-torus, in close agreement with bandstructure cal- 
culations. In Fig. 1 we present a direct comparison of 
the B C 2 data by Lyard et al [3] with our calculation for 
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Fig. 1. Upper critical field as a function of temperature for 
field in c-axis and ab-planc direction. Comparison with Lyard 
ct al. [3] (symbols). The inset shows the strong temperature 
dependence of the anisotropy ratio I ^c2- 

interband pairing strength r\ — 0.121 and c-axis disper- 
sion e c = 0.182 of the cylindrical Fermi surface sheet 
as described in Ref. [7] . For field in c-axis direction the 
temperature dependence follows very much a conven- 
tional one-band temperature dependence [8]. However, 
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Fig. 2. Zero energy density of states as a function of magnetic 
field in c-axis direction (solid line) and comparison with results 
from specific heat measurements by Bouquet ct al. [4] (symbols) 
on MgB2 single crystals taking B^ 2 — 4.2 T. The inset shows 
the field dependence of (A/T c ) 2 for the two gaps. 

a clear upward curvature close to T c is obtained for 
field in ab-plane direction. This results from a strong 
temperature dependence of the distortion of the vor- 
tex lattice, which at low temperatures is dominated by 
the (T-band, while close to T c both bands contribute 
equally [7]. 

Since our model offers a good description of _B C 2 in 
MgB2 we extended our calculations now to magnetic 
fields below B C 2 using the analytical method described 
in Ref. [9] for the calculation of the density of states 
within the vortex state. Within this method the field 
dependence of the partial density of states of each band 
at the Fermi level for field in c-axis direction is given by 
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Here, VF, a ,ab is the ab-plane component of the Fermi 
velocity of band a £ {a, n}. N a (0) is the normal state 
density of states, <E> U the flux quantum, and (■ ■ -) a de- 
notes a Fermi surface average over band a. For the 
cylindrical and torus Fermi surfaces considered here 
these averages can be done analytically. The field de- 
pendence of the gap A a (B) is obtained from a solution 
of the two by two gap equation in the vortex state [9]. 

In Fig. 2 we show the magnetic field dependence of 
the total density of states at the Fermi level (solid line) 
for the same set of parameters as in Fig. 1 showing fair 
agreement with the specific heat data of Bouquet et 
al. [4]. While the contribution from the large gap (a 
band, dotted line) to a good approximation follows the 
behavior of a single band s-wave superconductor, the 
contribution from the small gap band, dashed line) 
displays a rapid increase. The inset shows the field de- 
pendence of the square of the two gaps. It shows that 



the square of the large gap to a very good approxima- 
tion is a linear function of the magnetic field, while the 
small gap shows some reduction at higher fields. If we 
nevertheless approximate the squares of the gaps as 
linear functions in B/B^ 2 we find the following useful 
analytical expressions: 
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Here, E(k) and K(k) are the complete elliptical in- 
tegrals and k — 0.25 the ratio of the two radii of 
the torus Fermi surface as defined in Ref. [7]. At low 
magnetic fields from Eq. (1) we can obtain the fol- 
lowing limiting expressions: N a {B) ~ 1.67- 



and N n (B) w 1.01- 



These expressions 
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show that the rapid increase of the density of states of 
the 7r-band is due to two reasons: the small value of 
A-n (0) /T c and the higher in- plane Fermi velocity of the 
7r-band and thus has a different physical origin than 
the d-wave behavior known as Volovik's law [10]. We 
note, however, that these analytical expressions are 
only valid for field directed into c-axis direction. For 
other field directions the distortion of the vortex lattice 
has to be taken into account and we find important 
modifications. Details will be published elsewhere [11]. 
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